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Abstract: 

The statistical equilibrium properties of the linear sigma model are studied, with a view 
towards characterizing the field configurations employed as initial conditions for numeri- 
cal simulations of the formation of disoriented chiral condensates in high-energy nuclear 
collisions. The field is decomposed into its spatial average, the order parameter, and the 
fluctuations, the quasi-particles, and enclosed in a rectangular box with periodic boundary 
conditions. The quantized quasi-particle modes are described approximately by Klein- 
Gordon dispersion relations containing an effective mass that depends on both the temper- 
ature and the magnitude of the order parameter. The thermal fluctuations are instrumental 
in shaping the effective potential governing the order parameter, and the evolution of its 
statistical distribution with temperature is discussed, as is the behavior of the associated 
effective masses. As the system is cooled the field fluctuations subside, causing a smooth 
change from the high-temperature phase in which chiral symmetry is approximately re- 
stored towards the normal phase. Of practical interest is the fact that the equilibrium field 
configurations can be sampled in a simple manner, thus providing a convenient means for 
specifying the initial conditions in dynamical simulations of the non-equilibrium relaxation 
of the chiral field. The corresponding correlation function is briefly considered and used to 
calculate the spectral strength of radiated pions. Finally, by propagating samples of initial 
configurations by the exact equation of motion, it has been ascertained that the treatment 
is sufficiently accurate to be of practical utility. 
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1 Introduction 



The possibihty of producing and observing disoriented chiral condensates in high- 
energy colhsions of hadrons and nuclei has stimulated considerable interest over the 
past few years (for a recent review, see ref. 0). The basic premise is that the 
collision generates an extended domain of space within which chiral symmetry is 
approximately restored. If this happens, macroscopic pion fields may be generated as 
a consequence of the subsequent non-equilibrium relaxation towards the normal state. 
Such isospin-aligned domains may manifest themselves in anomalous pion emission 
0, 1^ of the type seen in the Centauro cosmic ray events 0. The experimental 
exploration of this phenomenon is of fundamental interest because it has a direct 
bearing on the mechanism of spontaneous chiral symmetry breaking. Efforts are well 
underway to search for the associated pion multiplicity fluctuations in proton-proton 
collisions 0. 

In order to assess the prospects for such a phenomenon to actually occur and be 
detectable above the background of other pion production processes, it is necessary to 
perform extensive dynamical calculations. This is a daunting task because the chiral 
degrees of freedom should be properly embedded in the complicated environment gen- 
erated in a high-energy collision which evolves from being primarily partonic at the 
early stage to entirely hadronic in the course of the chiral relaxation process. Fortu- 
nately, the study of suitable idealized scenarios can yield valuable insights regarding 
the prospects for observing the effect. 

Most dynamical studies have been carried out within the framework of the linear a 
model, treated in the mean-field approximation |^, |TD|, [T^, |T^, |T^, |T3|, [T^, [T^ . The 
physical problem divides naturally into three parts: 1) specification of the initial state 
in which chiral symmetry is approximately restored; 2) the dynamical evolution during 
which the chiral field cools and evolves towards the neighborhood of the physical 
vacuum; and 3) the analysis of the asymptotic state in terms of physical observables. 
The present work is particularly relevant for the first issue, as it provides a well- 
based framework for characterizing the chiral field configurations that are employed 
as initial states in the numerical simulations. In particular, a convenient method is 
developed for sampling the statistical equilibrium ensemble of chiral configurations 
and is expected to be of practical utility. 

After a brief reminder of the relevant features of the linear a model (sect. 0), we 
describe how thermal equilibrium can be treated approximately by means of a stan- 
dard linearization procedure (sect. ^. We then discuss and illustrate the statistical 
distribution of the average field strength, the order parameter (sect. and subse- 
quently turn to the properties of the quasi-particle degrees of freedom associated with 
the spatial field fluctuations (sect. An impression of the validity of the treatment 
is then gained by evolving samples of field configurations by the exact equation of 
motion (sect. ||). Finally, a concluding discussion is given (sect. ^. 
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2 The linear a model 



To set the framework for the subsequent developments, we start by briefly recalling 
the most relevant features of the formal framework. The present study is carried out 
within the mean-field approximation where the quantum field operators are replaced 
by their expectation values, thereby bringing the treatment to the level of classical 
field theory. This simplified treatment is expected to suffice for exploratory calcu- 
lations. Naturally, the mean-field approximation is only a first step towards a more 
complete description. 

The basic object of study is then the chiral field 

= (T(r,t) + zr-7r(r,t) , (1) 

where the three elements of the vector r are the Pauli matrices. Here the scalar field 
a{r,t) and the vector field 7r(r,t) are both real and can conveniently be combined 
into the 0(4) vector <j) = (a, tt). 

In the linear a model [|r^ , one introduces a simple local effective interaction energy 
density, 

1 



V{cf>) 



.2\2 



Ha 



(2) 



where (f) denotes the magnitude of the 0(4) vector (/> 

j=0 



a +77-77 



(3) 



with j = 0, 1,2,3 referring to the four chiral directions. The interaction V contains 
three parameters: A, the strength of the symmetric term; v, the location of its min- 
imum; and H, the strength of the symmetry-breaking term. As is commonly done, 
these parameters are fixed by specifying the pion decay constant, f-,, = 92 MeV, and 
the meson masses, = 138 MeV/c^ and rria- = 600 MeV/c^, leading to 
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20.14 
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) J IT 



86.71 MeV 



= (120.55 MeV)' 



(4) 

(5) 
(6) 



The precise values of the model parameters are immaterial, both in view of the simple 
nature of the model and, in particular, within the context of our present idealized 
study. The value for the hypothetical a mass is the most commonly adopted one; 
for m^r we have used the weighted average of the three observed pion masses, and 
for the pion decay constant we have simply taken two thirds of the pion mass which 
gives a simple value within the range of experimental data (the current values being 
= 92.4 ± 0.26 MeV and /^o = 84 ± 3 MeV [O]). 
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The Lagrangian density is given by 



C{r,t) 
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(7) 



The corresponding energy density is then 
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where the time derivative ip = hdt4> is the canonical conjugate of the field strength. 

The action generated by a given time evolution is given by iS = / drdtC{r,t). 
By demanding that S be stationary with respect to arbitrary variations of both the 
field strength 0(r, t) and its first derivatives, one obtains the associated equation of 
motion, 

□ </) + A(0^-t;^)0 = /7a- , (9) 

where □ = h^d^ — h'^c^A is the d'Alambert differential operator, and <t denotes the 
unit vector in the a direction. Since the equation of motion is of second order in 
time, the evolution of the system is fully determined once the initial values of the 
field strength 0(r) and the associated time derivative ip{r) have been specified. 



3 Approximate treatment 

We now discuss how the chiral system can be treated in a convenient approximate 
manner. For this purpose, we confine the system within a rectangular box and impose 
periodic boundary conditions. 

3.1 Linearization 

Our starting point is the equation of motion for the field strength. It is natural 
to decompose the field, 

0(r,t) = 0(t) + 50(r,t) , (10) 

where the first term is the spatial average, 4> =< (p >, so that < 6(p >= 0. The 
average part of the field, is often referred to as the order parameter, whereas the 
spatial fiuctuations, S(f){r), may be considered as quasi-particle degrees of freedom 



representing elementary excitations relative to the constant field configuration [|10 
The 0(4) direction of the order parameter plays a special role and it it convenient to 
employ a correspondingly aligned reference system, i.e. a system in which cf) = {(pQ, 0) 
and 64> = {6(f)\\, 6<p±). The thermal average -< 6(p S4> >- is then a diagonal 4x4 
tensor, as can be demonstrated by elementary means. 

Taking the spatial average of the full equation of motion (H) leads to an equation 



for time evolution of the order parameter ||20|| , 

h^d^fp + X[^l+ -< 64>^ 6^1 y -t;^]^ = Ha- , (11) 
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In deriving this equation, we have replaced the spatial averages < • > by thermal 
averages -< ■ as may be justified for a system of dimensions larger than the corre- 
lation length (the fiuctuations in different regions of the system are then independent 
and can therefore be assumed to reflect the thermal distribution from which the par- 
ticular fleld was sampled). We shall use -< >- to denote the diagonal elements of 
the isotropic 3x3 tensor < 6<p± S(f)± >, so that -< Scj)"^ >~=^ > — 1-3 ^ 6(j)]_ >- is 
the total fluctuation.^ Furthermore, < 50^ Sep > vanishes by symmetry. 

By subtracting the equation of motion for the order parameter, eq. (0), from the 
full equation of motion (|^), it is possible to obtain approximate equations for the fleld 
fluctuations, 

□ + /ijc^]50ll = 0, (12) 
□ + /iic^]50^ = 0, (13) 

where the effective masses /iy and fi± are determined by the auxiliary gap equations, 

/ijc^ = A (302+ ^ 502 ^ +2 ^ (50j5 >- -v^) , (14) 
fj^lc^ = A (4>l+ -< 502 ^ +2 ^ 501 y -v^) . (15) 



In arriving at this result, we have replaced products of two individual fluctuations 
by their respective thermal ensemble average values, 6(f)j6(f)j' ^ -< 6(f)j6(f)j' >~, and, 
furthermore, products containing three fluctuation factors have been contracted in 
the usual manner ||15|| , 



5(j)jS(j)j'6(f)j'i ^ -< 6(j)j6(f)j' >- 6(f)j" + ^ 6(f)j'6(f)j'i >- 6(j)j + -< 6(j)ji'6(j)j >- S(j)ji . (16) 



Equations ( P^ and ([TB|) describe the fluctuations as independent quasi-particles 
having the respective effective masses fi\\ and fi±, which in turn are given in terms 
of the fleld fluctuations, eqs. ([T^ ) and (p!5|) . This self-consistent relationship can be 
stated explicitly by invoking the expression for the associated thermal equilibrium 
fluctuations, 

h^C^^fl 1 T I /i||c2 

-= = IT? j^?^ " S5 ''"'^i; . (17) 

and similarly for -< 6(f)\ >~. The system is enclosed in a rectangular box with the 
volume Q = L^LyLz and periodic boundary conditions are imposed. The wave vec- 
tors k are then quantized {k^L^ = 2ttKx with = 0, ±1,..., etc.), so that the 
individual quasi-particle modes are enumerated by {K^, Ky, K^). The corresponding 
quasi-particle energy is determined by the dispersion relation el = tL^c^k"^ + fx^c^, 
where the effective masses /iy and n± lead to the energies e| and e^, respectively. Fi- 
nally, the prime on the summation sign indicates that the mode having k = should 
be omitted since it represents the order parameter rather than a spatial fluctuation. 

^We call attention to the fact that in the present notation -< (50j>- denotes the variance of the 
field fluctuations in one of the transverse directions, whereas in refs. flOj , |l4[ it represents three times 
that amount, namely the sum of the variances in all three perpendicular directions. 
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In eq. (p!7[), the symbol x indicates the thermodynamic hmit L ^ oo in which 
the quasi-particle modes form a continuum. The fluctuations can then be expressed 
analytically in terms of the Modified Bessel Function Ki, as indicated. Generally, 
the fluctuations decrease as the effective mass increases, for a fixed temperature T, 
and -< S(j)'j >~= when the effective mass fij vanishes. The global approximation 

-< 6(f)] (TV12)(1 + iic^/2'kT) exp(-/icVT) is good to better than 2%. 

The treatment above has been carried out under the simplifying assumption of 
thermal equilibrium which suffices for our present purposes. However, it is interesting 
to note that the results can be readily generalized to non-equilibrium scenarios by 
simply replacing the thermal averages -< ■ by the corresponding spatial averages 
< ■ >. In this manner, the system of equations would be closed and a conceptually 
simple dynamical description emerges. However, since the direct interaction between 
the quasi-particles is then ignored, any relaxation can occur only via the mean field, 
and the accuracy of such an approach should therefore be carefully assessed. 



3.1.1 Symmetries 

It is important to note that the equation of motion (pJ]) for the spatial average of the 
field strength, the order parameter 0, has four components, one for each of the chiral 
directions. 

In the special case when the symmetry-breaking term vanishes, H = 0, the in- 
teraction term in (0) depends only on the magnitude of the order parameter, 0o- 
The evolution of (f) is then equivalent to that of a particle in a central field in four 
dimensions and, accordingly, the dynamics separates into one radial and three an- 
gular modes. The rotational invariance implies that the angular motion is massless, 
a manifestation of the Goldstone Theorem (see, for example, ref. |^T|). The con- 
served currents associated with this 0(4) invariance are the Noether iso- vector current 
Vfj, = TT X dfj,7v and iso-axial-vector current = Tzd^a — ad^Tz (related to rotations 
around the a axis and rotations of the a axis, respectively). In the general case, when 
H does not vanish, the model remains invariant with respect to rotations around the 
a axis and the associated 0(3) symmetry ensures the conservation of isospin. 

We also note in this connection that the quasi-particles modes all have a finite 
effective mass (except right at the critical point). This should pose no puzzle since 
they have always finite momenta, whereas the above modes carry no momentum. 



3.2 Effective masses 



Utilizing the result (|T^, the coupled equations ([T^-[15D for the effective masses can 
be solved for specified values of the temperature T and the magnitude 00, provided 
that these parameters are sufficiently large, and we then have /iy > /i_L > 0. It should 
be noted that the effective masses are independent of the parameter H, since the 
symmetry-breaking term Ha is linear.^ 

^This simplifying feature holds only when A, v, and H are considered as the primary model 
parameters. When H vanishes it is often customary to readjust the other two parameters, A = 
m^c^ /2f'^ — 21.27 and v — f-^ — 92 MeV, and then the effective masses change correspondingly. 
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In fig. |l| we show the resulting effective masses and as functions of (po, 
for temperatures T up to well above critical. At any temperature, there is always 
a physical solution to the coupled equations ( p!4| - p!5| ) for the effective masses when 
00 > V. This is easy to see from eq. (^): At T = 0, when the fluctuations vanish, 
we have fj?j_ = A(0q — f^) and so //^ vanishes at 0o = "^^ and is positive for larger 0o; 
an increase of T will always increase the fluctuations, and hence the mass. Moreover, 
we always have /iy > fi±. 

Since the field fluctuations and the magnitude of the order parameter contribute 
to the effective masses in qualitatively similar ways, an increase of the temperature 
(and thus the fluctuations) will permit a further decrease of (po, so that the point 
at which vanishes is moved to ever smaller values of 0o- The behavior of fi\\, 
considered as a function of 0o, is nearly independent of temperature (as long as T is 
subcritical), except that each curve terminates at the point where the corresponding 
H± vanishes. This limiting curve is indicated by the dotted curve on the interval (0, v) 
and it is elementary to calculate, 

2 - - Ir^- ^ ^ ^ y^_Il(^5 + e-'^ii^'/'^) , (18) 



^0 



12 12 



/ijc^ = X{2v^ - T^) . (19) 
When the temperature reaches the critical value, T = Tc, there is a solution to 



eqs. ([14^151) for all values of 0o; and both effective masses vanish for 0o = 0. Then 



^ 502 and so = 2v'^, i.e. = 122.63 MeV with the adopted parameter 

values. The degeneracy in the masses, = /i_L, is a general consequence of the 0(4) 
rotational symmetry that emerges for 0o = and it therefore remains as T is further 
increased, with the common mass value /io increasing steadily. Since the effective 

mass at 00 = is given by fi^c^ = A(6 -< 50^ > f^), it becomes proportional to T 

at high temperatures, (Iqc^ ~ 1.59 T for T 3> v. 

The results displayed in fig. |I| have been calculated in the thermodynamic limit, 
L — s> oo, in which the quasi-particle modes form a continuum. However, the systems 
of interest in connection with high-energy nuclear collisions have finite volumes. The 
quasi-particle energies are then discrete and, as a result, the effective masses are 
modified. This effect is largest when the effective mass is small, i. e. near the critical 
conditions, because the absence of the constant mode in ( |17|) is then most significant, 
and it generally leads to smaller masses. Figure ^ shows the temperature dependence 
of the effective masses for either 0o = or for 0o = (the value associated with the 
physical vacuum), as obtained either in the continuum limit or for a cubic box with 
a side length L of either 8 fm or 5 fm. 

For 00 ~ and temperatures near critical, the effective masses are significantly 
reduced. As a consequence, the associated critical temperature is increased to 154 
MeV for L = 8 fm and to 172 MeV for L = 5 fm. However, for temperatures several 
tens of MeV above these thresholds, the effect is relatively small. The effect also 
disappears quickly as the order parameter is moved away from zero and by the time 
it reaches its vacuum value, 0vac = /tt = 92 MeV, there is hardly any effect for even 
the smallest box (L = 5 fm). 
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The effect of tlie finite size on tlie fluctuations in tlie field strength is illustrated 
in fig. ||. The upper curves show the typical magnitude of the thermal fluctuations, 
{-< 6(j)'^ >-)^^^, as functions of 0o over a range of temperatures, for the same three box 
sizes as in fig. ^, while the lower curves show the continuum result for the contribution 
arising from fluctuations directed along the order parameter (p, -< 50| y^^'^. The 
fluctuations are remarkably independent of the volume, except in the critical region 
T ^ Tc and (po ~ 0, as just discussed. For 0o = the 0(4) symmetry implies that the 
fluctuations are the same in all directions so that each of the lower curves starts out 
at half the value of the corresponding upper curve. For larger values of 0o the relative 
contribution from fluctuations along (p decreases progressively, since fj,\\ grows much 
larger than (see fig. p. 

3.3 Energy 

It is instructive to examine the total energy of the system. It is obtained by integrating 
the energy density (H) over the volume of the box, 

s{t) = Jdrn{r,t) = ^]< ^;iV + ^^V(v0)' + r > (20) 

= Q{Eo + E^^ + 6V) . (21) 

The total energy is conserved in time when the field is evolved with the equation 
of motion (P). The energy is composed of the three basic contributions exhibited in 
the second expression: one arising from the time dependence of the field strength, 
another associated with the spatial variation of the field, and a third resulting from 
the interaction. In the last expression, these contributions have been reorganized into 
three different terms which we will now discuss in turn. 

The first term is the bare energy density, i. e. the energy density that would arise 
if the field had no spatial variation, cp = (p and ip = ip, 



Ea 



^^0 + ^{(fl-v'^f - H(l)o cosxo 



Ko + Vo. (22) 



We have here introduced xo to denote the angle of disorientation, i.e. the angle 
between the four-dimensional order parameter (p and the a direction. The spatial 
average of a is then = (j) ■ & = 0ocosxo; and the corresponding mean magnitude 
of the pion field is ttq = 0osiiiXo- In the last expression Kq denotes the kinetic part 
of Eq and Vq is the bare potential energy density. 

The second term in (pi]) denotes the energy density associated with the quasi- 
particle degrees of freedom, in the approximation that these form a gas of independent 
particles, 

^qp = TT^T.l<i^9t6<j)jY + {hcV6<j),f + fiy6<p^^> (23) 
n c I 

= 7^EMlrfV + (^W + /.,V)|rfY). (24) 
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We have here invoked the general Fourier expansions of (pir) and il^{r), 



4^{r.t) = E0k(^)e''■^ (25) 

k 

i^{r,t) = E^kWe^'-^. (26) 

k 

The expansion coefficients are generally complex and are given simply in terms of 
spatial averages, 0^ =< 4>{r) exp(— zk ■ r) > and Vk =< ''/'(^) exp(— zk ■ r) >. We 
note that the components with k = represent the order parameter, </)q = (p and 
t/>Q = ip. Furthermore, the fact that the chiral field is real imposes the relations 
0k = 0-k and = tp_^. 

The third term in (^) corrects for the fact that the quasi-particles in fact interact 
via the fourth-order term in V , 

5V = \[<5<P^> - 2 < >2 - 4 < 5<P] >2] ^ -^<5<P^> . (27) 

The last result can be obtained by using the Gaussian approximation to express 
< Scf)^ > in terms of squares of quadratic terms. 

The energy of the system can be calculated as a function of the order parameter, 
assuming that the quasi-particle degrees of freedom are in thermal equilibrium. The 
resulting energy density is 

Ei±,f) = Eo{^oAo,Xo)+ ^^qp^(0o)+ ^5Vi(Po)y=Ko + VT{(po,Xo) . (28) 
Here the thermal equilibrium value of the quasi-particle energy density is given by 



34^.(^) + ^i^.(^) 
n"^ T n T 



(29) 



for each of the four chiral directions j, with r = T/fijC^. It can be expressed ana- 
lytically in the continuum limit, as indicated. Furthermore, the thermal equilibrium 
value of the correction term is quite well represented by its Gaussian approximation. 



SV ^ -^[^ (502 ^2 + 2 ^ 50j ^2 ^ Q ^ ^ (30) 

In the last expression in ( pHD the energy has been split into the kinetic energy 
associated with the time dependence of the order parameter, Kq = ipQ/2h^c^, and the 
remainder which can be regarded as an effective potential, 

Vt{(Po,Xo) = Vo+ ^E^^y + ^5Vy . (31) 

This effective potential has been plotted in fig. | for the particular situation when 
the order parameter (f) is directed along the cr axis, for a range of temperatures T. 
The directional dependence of Vr is simply given by the term H(f)ocosxo/^^c^- At 
zero temperature the effective potential can only be displayed for order parameters 
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whose magnitude 0o exceed v, since only then can the effective masses be calculated, 



as explained in sect. p.2| : For smaller values of 0o some of the low-energy quasi- 
particle modes are unstable and so have no well-defined thermal equilibrium. As the 
temperature is increased, the domain of stabihty extends ever lower towards 0o = 0. 
The limiting curve connecting the end points of all the subcritical curves is determined 
by the condition that vanish, implying /i| = X{2v'^ — T^). The result is shown on 
the figure as the dotted curve between (po = and 00 = "^^- 

For H = (fig. the effective potential has 0(4) symmetry and so the minima 
form the surface of a 4-sphere; its intersections with the a axis are indicated by the 
two arrows. The addition of the linear symmetry-breaking term (fig. ^) tilts the 
potential towards the a direction and, consequently, there is only one minimum for 
any given temperature and it is located on the a axis. 

For all values of H, the minima of the effective potential move gradually towards 
larger values of the order parameter. In this respect, the results are similar to the 
behavior of the energy of nuclear matter as a function of the density p (see, for 
example, ref. [^). However, while the energetics alone thus favors ever larger values 
of either 0o or p as the temperature is raised, the proper incorporation of the phase 
space (by means of the appropriate entropy) turns the picture around, as we shall 
now see. 



4 The order parameter 

We can now consider the distribution of the chiral field in statistical equilibrium at a 
given temperature T. Since the equation of motion (^ is of second order in time, we 
need to specify both the field strength <p{r) and its time derivative ■0(r) = hdt4>{r) 
in order to fully characterize a state. As we have already exploited above, the quasi- 
particle degrees of freedom can be regarded as a Bose-Einstein gas having thermal 
distributions determined by their respective effective masses, in the appropriate do- 
main of 00. Their behavior can then be calculated once the order parameter has been 
specified. 

4.1 The partition function 

It is natural to start the discussion of statistical properties by considering the partition 
function for the chiral field. It is obtained by summing the canonical weights for all 
the possible states. 



n 



J d^± e-^(^o+^^^^)Zqp = j c/V Wt{±,±) , (33) 



where £ is the total energy of the system given in eqs. (|20| - piD . The functional integral 
over all states {il)[r), (f)[r)) has first been expressed as a multiple integral over the 
four- dimensional Fourier amplitudes -j/j^ 0k (^^^ ^l^- (llHl))- If ^^e integration 
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over the field fluctuations {i.e. those amplitudes having k 7^ 0) is performed, the 
remaining integrand is the statistical weight Wt expressing the relative probability 
for finding the system with a specified value of the order parameter. In order to 
derive an expression for Wt, we have replaced the correction term 6V (eq. (p7D) by 
its approximate thermal average -< SV >-, for each value (i/j, 0) of the order parameter 
(it depends only on 0o)- The integral over the quasi-particle degrees of freedom can 
then be performed and yields the quasi-particle partition function Z^p. 

This quantity factorizes into contributions from each of the four chiral directions, 
2qp = where 



n' 

k 



k 



27r 



(34) 



Since the complex Fourier amplitudes and (f)^' are subject to symmetry relations 
under k — > — k, as already noted above, is it more convenient to use the representation 
in terms of trigonometric functions, eqs. (^) and (^). The different wave numbers 
are then fully decoupled, and for each k the integration variables are the (positive) 
amplitude Ck and the associated phase ?7k, as indicated in the last expression above. 
The phases are unimportant here since they each merely contribute a factor of 27i. 
Moreover, the combination rik = ekC^/fi^c^ can be interpreted as the average number 
of quanta in the mode k (see eq. (0)). Since we want to take account of the quantal 
nature of the quasi-particle degrees of freedom (see |5.1.1| ), we replace the continuous 
integral over rik by a discrete sum over the possible integer values of the occupation 
number, thereby obtaining the standard expression for an ideal Bose-Einstein gas, 



U) 



qp 



n 

k 



nk=0 



k k 



(35) 



We have used the fact that the quasi-particle energy is additive, = I]k"-kefe- 
Furthermore, = 1 + where fk = (e'^'^Z^ — 1)^^ =-< rik >- is the mean number 
of quanta (in the chiral direction j) having a wave number of magnitude k, and the 
corresponding energy is determined by = h^c^k'^ + jJ^'^jC^- 

It follows that the free energy density of the quasi-particles is Fqp = Y,j F^p with 



qp 



--lnZ(j) 



E' [^kfk - T{fk In fk - fk In fk 



^ qp ^ 



-Ts: 



(36) 



In the resulting expression, -< E^J >- is the average thermal energy density of the 



quasi-particle modes associated with the chiral direction j, eq. 



), and s!^^ is their 

entropy density. The total entropy density at the temperature T is then S*!- = Sf. + 
bf)) represents the entropy density of the field fluctuations along the 



35*^, where Sj.{ 
direction of (/> p3 



4 = ^E' [4' In 4' -/I In /i 



de 

:2 2^ 



e e- 



/ijc^)2 



/(e) In /(e) -/(e) In /(e) 



(37) 
(38) 
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and the contribution to the entropy density from fluctuations in a perpendicular 
direction, Sj^{(j)Q), can be calculated analogously by replacing /iy with 

It is evident from fig. ^ that the energetics favors rather large magnitudes of the 
order parameter, as would be expected from the form of the interaction potential V. 
However, the density of quasi-particle states is higher at small values of 00, since the 
effective masses are then smaller, for a given temperature.^ This interplay between 
energy and phase-space is quantified in the statistical weight (see eq. (P^)), 

Wt{±,^) = J V[iP', (/)'] e-^[^''^']/^ - V^) - ±) (39) 

We have used the approximate decomposition (^) of the total energy which causes 
the partition function to factorize and, in the last expression, we have introduced the 
free energy density for the order parameter, Ft(0) = Vt(0o5 Xo) ~ 7'S't(0o)- 

The free energy density is illustrated in fig. |^ for a range of temperatures above 
which the trend is obvious. Since the entropy density of a free gas at fixed T goes 
up when the particle mass is reduced, the entropy provides a restoring force towards 
symmetry (recall that a reduction of the order parameter 0o leads to smaller effective 
masses). As a result, for H = and temperatures well above Tc the free energy 
density then has its minimum at 0o = 0. As the temperature drops, the symmetric 
minimum grows ever more shallow. Near T ^ 175 MeV a secondary minimum appears 
at 00 ~ 50 MeV, and it becomes the lowest one at T ?a 171 MeV. As T decreases 
further the minimum gently approaches (pQ = v . The abrupt change in the location of 
the free-energy minimum is characteristic of a first-order phase transition. In contrast, 
when H takes on a realistic finite value (see eq. (j^)) the minimum of Ft is always 
located on the positive part of the a axis and it moves smoothly outwards as the 
system is cooled, with no drastic change at any temperature. The resulting behavior 
is shown by the solid curve in fig. ^ 

4.2 Distributions 

The preceding analysis implies that the statistical weight Wri'tp, <p), which is defined 
on the eight-dimensional phase space of the order parameter, in fact depends on only 
three quantities: the speed ipo = \'^\, the magnitude <^o = 101, and the disalignment 
angle xo- Moreover, since the entropy depends only on the effective masses, which in 
turn are determined solely by (po, the dependence of the free energy on the direction 
of arises only through the symmetry-breaking H term in the interaction V (eq. 
(^). Consequently, the corresponding projected probability distribution factorizes, 

Pt{A,(Po,Xo) = /rfV/rfV5(^^-^o)5(0o-0o)5(xo-Xo)W^T(:0'>') 

= P^(V^o) P<p{<Po) Pxi^o;Xo) ~ ^0 00 sin^Xo W^t(^o, 0o, Xo) • (41) 

•^The importance of the quasi-particle entropy in determining the statistical properties of the 
model was already emphasized by Bardeen and Moshe a decade ago. 
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Here ip^ and 0q are the Jacobians associated with the transformation from the four- 
dimensional vectors t/> and (f) to their magnitudes ipQ and respectively, and sin^ xo 
is the Jacobian arising from the three-dimensional nature of the space perpendicular 
to the order parameter (f). Furthermore, the three normalized probability densities 
are given by 

P^^o) = ^0 ^^P[-^^] ' (42) 

P0(0o) = -A/; 0^ exp[-^](^V^T(0o, Xo = 0) - 5t(0o))] , (43) 

Vt H 

Pxi(Po;Xo) = AT^l^o) sin^xo exp[-^3-j— 0o(l - cosxo)] • (44) 

For the normalization constants we have A/"^^ = 2k'^, with the convenient abbreviation 
K = Th^cyn, and = (vr/0/i(0 exp(-0, with ^ = H(j)o/K. 

The kinetic part, P^{iPq), is the projection of an isotropic four- dimensional normal 
distribution having a total variance oi -< ipQ >~= 4k. The average kinetic energy 
associated with the time evolution of the order parameter is then -< K >~= 4T, as 
would be expected in four dimensions. The mean speed is ipo >~= y^Ovr/t/S. 

Considered as a function of the four-dimensional order parameter 0, the probabilty 
density has a maximum where the free energy has a minimum, as is evident from eq. 
(pQl). However, the jacobian factor ~ 0q associated with the projection from (f) to 
00 increases the most probable value of the magnitude 0o- This effect depends on 
volume Q, as is illustrated in fig. |^. The results for a finite box with L = 8 fm (dashed 
curve) are nearly identical to those obtained for L ^ oo, while L = 5 fm (dots) leads 
to a more significant increases in 0o at the higher temperatures. The most rapid 
increase of 0o occurs near T ^ 200 MeV which is well above the critical temperature 
Tc (arrow) at which the effective masses drop to zero. 

The directional distribution Px{<p0i Xo), which depends parametrically on the mag- 
nitude 00; grows broader at high temperature where 0o decreases. In the temperature 
range considered here, up to several hundred MeV, Pxixo) is well approximated by 
the projection of a three-dimensional normal distribution with a total variance of 3/^, 

and the average magnitude of the disalignment angle Xo is then ^ Xo y's/vr^. 

Although the locations of the minima in the free energy appear qualitatively dif- 
ferent when comparing H = with H > (figs. ^ and respectively), the corre- 
sponding curves for the most probable 0o are rather similar when finite volumes of 
nuclear size are considered, due to statistical fluctuations. This is also illustrated in 
fig. 1^. The thin curve is analogous to the solid curve and shows the dependence of 0o 
in the thermodynamic limit. (Recall that the vacuum value is 0o = f = 86.71 MeV 
for H = 0.) At T f» 171 MeV the value changes abruptly. The open dots are ana- 
loguos to the solid dots and show the most probable 0o for the box having L = 5 fm. 
The bias by the jacobian factor 0g (see eq. (^3D) causes the behavior to become quite 
similar to that obtained for H > 0. In particular, there is no discontinuity near 
T 171 MeV where the minimum in the free energy jumps between 0o = and 
f» 56 MeV. We finally note that for relatively low temperatures the Jacobian causes 
the most probable value of 0o to slightly exceed the vacuum value. 
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4.2.1 Distribution of the order parameter 

For any finite volume Q, the order parameter the exhibits fluctuations around its 
most probable value. The full width at half maximum of the 0o distribution is shown 
by the horizontal bars in flg. |^ for the smallest box (L = 5 fm) where the fluctuations 
are the largest. The relative smallness of the fluctuations indicates that the order 
parameter is distributed within a rather limited range of values. This feature is 
further illustrated in flg. ^ which depicts the entire distribution, P^i^cpo), over a range 
of temperatures, for a cubic box with side length L = 8 fm, which is our standard 
scenario. The solid curves in flg. |^ have been obtained by scaling the continuum 
results to the flnite volume. The dashed curves show the effect of properly quantizing 
the problem {i.e. summing over the discrete modes rather than integrating). The 
effect is very small, because the distributions are peaked well inside the respective 
domains of stability. 

For small temperatures, the distribution is narrowly peaked near the vacuum value 
0vac=/7r (indicated by the arrow). As the temperature is increased, the distribution 
broadens and gradually begins to move inwards towards smaller values of 0o- The 
width of the distribution increases from zero at T = 0, exhibits a maximum near 
T 220 MeV, and then slowly shrinks as T ^ oo. The somewhat counterintuitive 
decrease of the fluctuations at high temperatures is due to by the fact that the growing 
thermal fluctuations cause the interaction to become progressively more repulsive, 
causing the effective potential Vt to become ever more conflning. 

We recall that for subcritical temperatures, T < T^, our treatment can only be 
carried out above a temperature-dependent minimum value of 0o- As it turns it, this 
principal limitation is of little practical import. For example, for T = 80 MeV (the 
lowest temperature shown in flg. ^ the distance from the centroid (at ~ 91 MeV) 
to the boundary (at ~ 69 MeV) is over ten times the dispersion of the distribution 
[k, 2 MeV). Thus, for any temperature, the distribution P^[(j)Q) is sufficiently narrow 
to make incursions into the respective unstable regime extremely unlikely. 

A more global impression of the statistical distribution of the order parameter 
(f) can be obtained by considering a contour diagram of the projected probability 
density P(0O! Xo) = -P</>(0o) Pxi'^'o', Xo) ■ Such a plot is displayed in flg. | for a box 
with L = 8 fm. The abscissa is the projection of the order parameter onto the a 
axis, (Jo = ■ 6" = 00 cos xo 5 and the ordinate is the magnitude of its perpendicular 
component, ttq = 0osinxo- For each temperature, the dots show the location of the 
maximum and the half- maximum contours are traced out. The solid contours and 
dots refer to the continuum treatment, whereas the dashed contours and the open dots 
(shown for two temperatures only) are obtained with a quantized treatment, which 
is seen to have little effect. The proflles in flg. |^ are the projections of P(0o, Xo) onto 
the magnitude (p^. It is evident from the above results that the 0(4) symmetry is 
far from restored at the temperatures expected in the planned high-energy nuclear 
collisions. Instead the order parameter is in fact distributed within a fairly limited 
domain, with both its magnitude 0o and its angle of disalignment xo exploring only 
rather narrow ranges. 
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4.2.2 Distribution of the effective masses 



The effective masses are functions of the magnitude of the order parameter, fi\\{(f)o) 
and /i±(0o)- Consequently, the thermal fluctuations in the order parameter will cause 
the effective masses to fluctuate as well. Since their distributions may be of some 
interest (see sect. ^.31) , we consider briefly the corresponding probability densities for 
the effective masses, 

P^(/ic2) = J rfV 5(/i±(0o)c' - /ic^) , (45) 

and analogously for P\\{fic^), where P{4>) = P0(0o)-Px('^o; Xo) is the probability density 
for the order parameter 0. In order to gain an impression of these distributions, we 
show in fig. |^ the most probable effective masses as functions of the temperature. 

The transverse effective mass increases steadily from its free value as the 
temperature is raised, whereas the parallel effective mass /iy drops by nearly fifty per 
cent before finally turning upwards. For temperatures above T ^ 300 MeV the two 
effective masses are practically degenerate and gently approach their asymptotic form 
fic^ ^ 1.59 T. We note that the effective masses always exceed the temperature, so 
it is never reasonable to ignore the effective masses. We also note that the transverse 
mass has practically no strength near the free pion mass, until the temperature is 
below Tc (see fig. p!3|) . 

The parallel effective mass yUy exhibits a pronounced minimum near T ^ 240 MeV. 
It is in this region that the order parameter exhibits its most rapid evolution with T 
(see fig. ^ and the temperature at which the minimum in /iy occurs may be interpreted 



as the effective critical temperature of the model |25 . 



4.3 Sampling of the order parameter 

The above analysis provides a convenient basis for sampling the order parameter in 
accordance with the statistical weight Wt{<P,iI^) given in eq. (|39D , and we describe 
briefly how this can be accomplished in a manner that is both quick {i.e. the compu- 
tational effort is small) and efficient {i.e. no effort is expended on rejection). 

As noted already, the time derivative ip is governed by a four- dimensional nor- 
mal distribution which is isotropic and entirely decoupled from the other degrees of 
freedom (see eq. (^)). Thus it is elementary to sample this quantity. 

The most complicated sampling concerns the magnitude (po, due to the intricate 
structure of its probability distribution, as discussed above. However, the numerical 
effort required is quite modest. The most efficient method requires a precalculation 
of the effective masses as functions of (po, for the particular T of interest. This is 
quickly done by proceeding as described in sect. [X^. The xo-independent part of the 
effective potential, Vt(0o!O), can then be obtained together with the corresponding 
entropy 5*^(00) • In effect, the probability distribution for 0o can be pretabulated, 
P<l>{4>o) (see eq. (|43|)), and it is then a numerically trivial task to sample 0o- 

Once the magnitude 0o has been selected, it is straightforward to sample the dis- 
alignment angle Xoy using either the exact form (^ij) or its Gaussian approximation. 
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In order to orient (f) in the tt subspace, there remains the task of selecting the re- 
maining 0(3) spherical angles "i^o and ipo, upon which the order parameter is given 
by = (00 cos xo, 00 sin Xo sin 'do cos v^o, 0o sin Xo sin 'do sin yjo, 0o sin Xo cos 'do) ■ 



5 The quasi-particles 

We turn now to the discussion of the quasi-particle degrees of freedom associated 
with the spatial variations of the chiral field, 6(p{r). Once the magnitude of the order 
parameter, 0o, has been chosen (see sect. |473| ), the quasi-particle degrees of freedom 
are fully characterized, via the effective masses /i||(0o) and /i±(0o), and it is then 
possible to sample them them appropriately. 



5.1 Sampling 

On the basis of the above developments, it is possible to devise a simple and efficient 
method for performing the statistical sampling of the chiral field, thus putting the 
initialization of dynamical simulations on a formally sound basis. 

Since the different quasi-particle modes can be regarded as effectively decoupled, 
the sampling is best done by making an expansion into the elementary modes, 

^011 (r, t) = ' cos(k -r-Jlt- 4) , (46) 

and similarly for the three transverse chiral components d4>j_{r,t). Here the energy 
ek = hcuk is determined by the Klein-Gordon dispersion relation, = fi^c^k"^ + jj^c^ 
(with /i being the appropriate effective mass, /iy or The phase ?7k is random in 
the interval (0, 27r) and is thus trivial to sample. Furthermore, the real (and positive) 
amplitude Ck can be related to the number of quanta nk by considering the energy 
carried by the mode. 



n 



(47) 



We have here omitted the zero-point energy [^], thereby eliminating the associated 
ultraviolet divergence. Although not quite correct, this approach is justified a pos- 
teriori by the apparent good quality of the resulting approximate treatment. (The 
dynamical tests discussed in sect. ^ are here very important, since the statistical 
samples would not remain stationary under time evolution if the treatment were sub- 
stantially wrong.) 

Thus the problem has been reduced to sampling the number of quanta n^. Since 
the probability for finding a particular number of quanta in the mode k is given by 
Pi^k) = (1 ~ exp(— efc/T)) exp(— nke^/T), it is elementary to sample the integer nk 
appropriately.^ We note that the thermal average of nk is equal to the occupancy fk 
employed in the calculation of the entropy, -< y= fk- 

"^The corresponding algorithm for this task is very simple because nt can be regarded as counting 
the number of successive times the sampling of a standard random number (i.e. uniform on (0, 1)) 
yields a value below exp(— efc/T). 
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Once the amplitudes and phases have been selected, the expansion P6|) readily 
yields the initial value of the field fiuctuations, 5</>(r, 0). The corresponding conjugate 
momentum, ^t/j = hdtScf), readily follows since (^61) implies 



^V^ll (r, t) = [^y Y^dcl sin(k -r-ult- 4) . (48) 

The entire state of the chiral field, {(f){r),il){r)), has then been selected at the time 
t = 0. When the equations of motion are propagated by a leap-frog method, the field 
strength S(f){r) is calculated at the times t„ = nAt while the momentum 6ip{r) is 
obtained at the intermediate times. The appropriate initial Sif){r) can then easily be 



obtained by evaluating the expansion (^8[) at t = At/2, after Ck and ?7k have been 
selected. 

Finally, since the sampling has been done in a system aligned with the order pa- 
rameter (f) (in which the mass tensor is diagonal), an 0(4) rotation of the sampled 
field configuration is required in order to express the state with respect to the chi- 
ral directions. This is readily accomplished on the basis of the angles (xo? ^o? V^o) 
specifying the direction of 0. 



5.1.1 Occupation numbers 

Since the quasi-particles represent bosonic modes, it is useful to know how large the 
occupancies can become. The thermal occupation numbers are given by -< rzk y= 
fk = l/(exp(efc/T) — 1), and since < fi\\ in any given scenario the largest occu- 
pancies occur for the transverse modes, > /|. Moreover, for a given temperature 
T, fk is largest when the momentum vanishes, /o = l/(exp (/ic^/T) — 1), which then 
provides an upper bound on fk (recall that /c > for the quasi-particles modes). 

Considered as functions of T, the bounds /g and f(^ start out from zero, display 
maxima well above T^, and then drop off towards a common constant value at high 
temperatures. The (common) limiting occupancy is /o ~ 0.26 since /zc^ ~ 1.59T when 
T ^ oo. (This feature is a direct consequence of the repulsive self-interaction of the 
chiral field and is in marked contrast to the ever increasing occupancy characteristic 
of free bosons in a thermal bath.) The maximum values attained by /g and Z^*" 
are about 0.38 and 0.48 and occur at approximately a temperature of 265 and 235 
MeV, respectively. Since these values are well below unity, the system is never very 
degenerate. 

Nevertheless, there is an important advantage to using Bose-Einstein rather than 
Boltzmann statistics for the field fiuctuations, as we will now discuss. When classical 
statistics is used, the occupation probability is Uk ~ T/ek when ek ^ T and there- 
fore the total quasi-particle density grows as the square of the upper limit on e^, a 
manifestation of the Rayleigh- Jeans divergence. Such a description would be entirely 
wrong in the present context and, moreover, it would be numerically ill-behaved. By 
contrast, the quantal occupation probability falls off exponentially and the density of 
quanta is finite. By adopting the Bose-Einstein treatment, we then get a much more 
realistic formal description of the statistical properties of the system and, in addition. 



16 



the numerical treatment becomes straightforward. Of course, if a configuration sam- 
pled on the basis of the quantal statistics is propagated for a sufficient length of time, 
it will eventually exhibit classical features, since the equation of motion is entirely 
classical. It is therefore fortunate that our treatment is expected to be applied only 
to processes that are far faster than the time scales assoicated with the reversion to 
classical statistics. In particular, for the formation of disoriented chiral condensates 
in high-energy collisions the relevant time scales are of the order of a few fm/c, while 
our numerical studies exhibit no ultraviolet run-away for at least several tens of fm/c. 



5.2 Correlation function 

It is particularly interesting to calculate the correlation function of the chiral field 
since this quantity determines the spectral distribution of the emitted field quanta. 

The density matrix associated with the quasi-particle degrees of freedom is a 4 x 4 
tensor, 

C(ri2,ti2) = -< 50(ri,ti) 5</)(r2,t2) >- • (49) 

The average is over the ensemble of systems considered, in the present case a thermal 
ensemble held at the temperature T. Since a system in equilibrium is invariant in time, 
the correlation function depends only on the time difference tu = ^1—^2- Analogously, 
the translational symmetry of the scenario implies that the spatial dependence is via 
the separation ri2 = Vi — r2- Moreover, to the extent that there is invariance under 
spatial rotations, only the magnitude ri2 = 1^12! enters. In principle, these spatial 
symmetries are broken when a finite box is considered but the effect is insignificant 
and can be disregarded in the present study. 

Utilizing the expansion (|46|), it is elementary to show that the correlation tensor 
C is diagonal with the elements C\\ and C_l, where 

^■^c^^/ 1 cos(k ■ r — cufct) 1 he f°° J sinfcr 

k 

with e| = h'^ujl = h^k'^c^ + /iyc"^. An analogous expression holds for C±{r,t). Since 
^^c^fk/^k is equal to the thermal average -< we recognize that the result (pOD 



^, , ft c%-^/ i cos k ■ r — cufci L he f°° , smfcr 

C|| r,t = ^— y ^ X / de ^ cosw^t , 50 



corresponds to the expression (146.10) given in ref. [27 
The usual correlation function is the trace of C, 

C =-<5ct){riM)-5(t>{r2M)^= tr C = C||(ri2, ti2) + 3Cx(ri2, ^12) • (51) 

Its overall magnitude is set by its value at zero, which is simply the corresponding 
fiuctuation in the field strength, C(0, 0) =-< 50^ >~. It is therefore convenient to divide 
C by this quantity and so define the reduced function C12 = C'(ri2,ti2)/ -< 
Since C12 is unity when ri = r2 and ti = ^2, it expresses the space-time attenuation of 
the correlation between the field-strength fiuctuations at different space-time points. 
In general, C12 ~ (l/ri2) exp(— /icri2//i) in the limit of large separations, ri2 00, 
so that /i//ic provides a simple measure of the correlation length. In the special case 
when the effective mass vanishes, the reduced equal-time correlation function is given 
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on analytical form, C12 x (3/C)(coth(^ — l/C)? where C, = TcTr 12 /he. In this extreme 
case, the correlation function falls off only as ~ 1/^12. 

Figure |10| shows the reduced correlation function Cu obtained at the critical tem- 
perature, T = Tc. Its appearance depends on the magnitude of the order parameter, 
00, through the effective masses. As (po is increased from zero to its vacuum value 
/,r = 92 MeV, the corresponding effective masses increase from zero to their free val- 
ues and the correlation function falls off ever more rapidly. The correlation function 
thus exhibits a significant sensitivity to the order parameter. We note in particular 
that for the most probable value, 0o ~ 87 MeV, the attenuation of Cu is consider- 
ably faster than for 0o = and its width has dropped by over a factor of two. It is 
therefore important to take proper account of the order parameter when calculating 
the correlation function. 

Figure |Tl| shows how the reduced correlation function evolves with temperature 
when the most probable magnitude of the order parameter is employed. At high tem- 
peratures the effective masses grow nearly in proportion to T so then the correlation 
length tends to zero, as is borne out by the steady shrinking of Cu. For temperatures 
below critical, the field fluctuations are predominantly associated with the transverse 
modes, since those have the smallest effective mass, fi± <^ and the correlation 
length grows ever larger. 

It is common to characterize the system by the "correlation length" , defined as the 
first moment of the equal-time correlation function which diverges when /x — * 0. For 
the present discussion, it is more convenient to characterize Cu by its full width at half 
maximum for equal times, ri2, since this quantity is always easy to extract, even when 
the mass vanishes. Figure |1^ shows this measure of the correlation length as a function 
of temperature, using the most probable value of 0o- For temperatures near and below 
T ^ 200 MeV the dominant fluctuations are perpendicular to the order parameter 
since the corresponding effective mass is relatively small. For higher temperatures 
the asymptotic regime is approached where the chiral symmetry is approximately 
restored and the fluctuations are similar in all four chiral directions. 

It is sometimes of interest to also consider correlation functions involving 6ip, the 



time derivative of the local field strength [|T^]. Those can be obtained in a similar 
manner. 



5.3 Radiation spectra 



The evolving chiral field may give rise to mesonic radiation, in analogy with the emis- 
sion of photons by a time dependent electromagnetic field. The rate for production 
of a field quantum having energy E and momentum p is proportional to the square 
of the corresponding Fourier amplitude ^ 

MN^ If. r 6(j)j{r, t) e 



E- 



dp 



dr 



\{p-r-Et) 



(52) 



Here j denotes the particular 0(4) direction considered so that j = gives rise to 
isoscalar a-like mesons and j > represent three components of the isovector pion-like 
mesons (with j = 3 corresponding to ttq, say). 
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A uniform system in equilibrium has both temporal and spatial invariance and 
the specific radiation rate, z/j(p, E) {i.e. the production per unit volume and per unit 
time), is then given by the Fourier transform of the correlation function 

1 rIN- r r 



where Q is the volume of the system and to is the time interval considered. The 
Fourier transform of the correlation function is easy to obtain in the continuum limit. 



C,(p, E) X h'c'^^^ S{m'c' - fiy) , (54) 

where we have assumed that the energy and momentum of the radiated mesons are 
related by E'^ = p^c^ + im?c^. 

The above result holds for a specified value of the order parameter 0o which 
determines the effective mass ^j. As we have discussed, the order parameter has in 
general a statistical distribution, P{4>o), giving rise to corresponding distributions of 



the effective masses, P\\{n) and P±{ij) (see sect. |4.2.2| ). The resulting production rate 
can then be obtained by integrating over the appropriate thermal mass distribution. 
Thus, for example, the specific production rate for ttq mesons is 



and the rates for the charged pions are similar .0 This result is easy to interpret: The 
spectral distribution of the radiated mesons is of thermal Bose-Einstein form, with 
the temperature given by the value characterizing the source itself, and the overall 
normalization of the radiation rate is proportional to the probability that the space 
and time evolution of the field matches the particular dispersion relation for the type 
of meson considered. 

As an illustration, fig. |1^ displays the strength function Pj_(m^c^) obtained by 
evaluating eq. (|^) for = m^r as a function of the source temperature. As the 
temperature is decreased towards zero, the centroid of the /i^ distribution moves 
down towards the free value rriTr and at the same time its width keeps shrinking 
towards zero (see fig. §). Through the supercritical regime the pion strength then 
grows approximately exponentially. Then a plateau is reached where the increase 
caused by the approach of the centroid to m^r is counterbalanced by the decrease due 
to the shrinkage. As a result, the strength is nearly constant from 90 to 30 MeV. 
Finally, after the centroid has practically reached m^, the free strength exhibits a 
rapid rise as the fi± distribution approaches a 6 function. On the basis of this result, 
one would expect radiation of free pions to be unimportant until the temperature has 
dropped below Tc. 

^The rates considered here pertain to the idealized scenario of a macroscopically uniform system 
(possibly enclosed in a periodic box) , and hence they differ from those describing the emission from 
a finite source into the surrounding vacuum. In particular, the familiar kinematic enhancement of 
the faster-moving ejectiles is absent. 
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6 Dynamical tests 



Once the initial state has been prepared, for example by means of the statistical 
sampling described above, the chiral field may be propagated in time by means of 
the equation of motion (^, which is straightforward to implement in either x or k 
representation. 

It is possible to exploit the dynamical evolution to test the validity of our approx- 
imate statistical treatment. If the system is ergodic, as would be expected because 
of its non-linear interaction, a dynamical trajectory will explore the space of possible 
field configurations in accordance with the appropriate microcanonical weight. Con- 
versely, an ensemble of field configurations that has been sampled statistically should 
not exhibit any change under time evolution. These features provide a convenient 
means for checking our treatment and we give two illustrations below 



6.1 Average field strength 

Perhaps most vividly, we show in fig. the early trajectories of the order parameter 
(f) for a sample of eight configurations, considering a box with L = 8 fm and either 
T = 200 MeV or T = 240 MeV. The dashed contours are those already given in 
fig. ^ indicating where the projected probability density has fallen to half its maxi- 
mum value and the centroids are indicated as well (open dots). For each individual 
trajectory, the initial location is indicated by the solid dot and the attached solid 
curve shows the trajectory up to t = 1 fm/c. The fact that the initial points reflect 
the calculated statistical distribution provides an elementary test of the numerical 
sampling algorithms. Less trivial is the fact that the dynamical trajectories indeed 
appear to explore the region predicted by the approximate statistical distribution. 

This correspondence can be made more quantitative by studying how the distribu- 
tion of the order parameter evolves in the course of time. This analysis is illustrated 
in fig. |T5|. The approximate distribution P^{(po) given in (|43|) is indicated by the 



solid curves for T = 180 and T = 240 MeV. Forty individual systems have then been 
prepared by sampling their field configurations as described above and they have sub- 
sequently been propagated by the equation of motion (H) up to the time t = 10 fm/c. 
In the course of the evolution, the value of (po is extracted at regular intervals and 
binned into slots that are 5 MeV wide. In this manner the time-averaged distribu- 
tion of 00 can be determined and the dashed curves display the result (which is not 
sensitive to an increase of either the maximum time or the sample size). The overall 
agreement with the initial distribution is very good. There is generally is slight shift 
outwards, amounting typically to 1-2 MeV, which suggests that our approximate ther- 
mal distributions are centered at correspondingly too low values of (pQ. Consequently, 
for applications where such a relatively minor displacement is unimportant, the test 
demonstrates the validity of our approximate method for obtaining the statistical 
properties of the model. 



^It is important to recognize that whereas the statistical properties have been obtained by assum- 
ing that the quasi-particle degrees of freedom are effectively decoupled, no such assumption is being 
made in the dynamics, since the trajectories are obtained by solving the full equation of motion (p|). 
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6.2 Field fluctuations 



Figure ^ displays the correlation function for the pion components, Cf2(si2), for a 
box with L = 5 fm prepared at T = 240 MeV. The short-dashed curves show the 
correlation function obtained on the basis of ten sampled configurations, and the 
long-dashed curves then indicate the corresponding result after those systems have 
been evolved up to the time t = 10 fm/c. The two curves that go up again have 
been obtained by aligning the relative separation along one of the cartesian axes, 
while it is directed diagonally for the other two curves (the periodicity is then ^/SL 
and so their eventual rise is only barely visible). For reference is shown the exact 
thermal correlation function for either the finite box considered (solid curves) or the 
continuum limit (dots). While this latter curve tends to zero for large separations 
(and in fact falls off monotonically) , the correlation function for a finite box drops to 
a negative value, because its spatial average must vanish. 

The correlation function remains remarkably invariant in the course of time. This 
indicates that our treatment, including the sampling procedure, in fact yields a good 
approximation to the correlated field fluctuations. As a quantitative measure, one 
may consider the full width at half maximum ri2. The continuum value is 1.33 fm, 
slightly larger than the thermal result for the finite box, 1.30. For the sample of ten 
initial configurations we find 1.256 and 1.271 for the cartesian and diagonal directions, 
respectively, which have evolved into 1.263 and 1.290 at t = 10 fm/c. So there is no 
significant change in the width of the correlation function. 



7 Concluding remarks 

The present work was motivated by the current interest in disoriented chiral conden- 
sates, particularly the various dynamical simulations carried out with the linear a 
model 0, ^ |Tl], |T3|, [1^, . Those calculations follow the non-equilibrium evolution 
of the cooling chiral field in order to ascertain the degree to which coherent domains 
develop. Since the dynamics is inherently unstable, with the low-momentum modes 
experiencing rapid amplification, one may expect a significant sensitivity of the re- 
sults to the initial conditions, with a commensurate degree of difficulty regarding their 
interpretation. Consequently, caution is required when characterizing the ensemble 
of initial field configurations employed. 

In order to provide a useful framework for this aspect of the problem, we have 
explored the statistical properties of the linear a model, by confining the system to a 
rectangular box held at a fixed temperature. Although this problem can be treated 
exactly ||31[ , we have found it preferable to linearize the equations of motion since our 



view is towards practical calculations. The resulting treatment then becomes very 
simple and appears to be sufficiently accurate in the present context. 

The problem separates into one concerning the spatial average of the field, the 
order parameter, and another dealing with the field fluctuations, referred to as the 
quasi-particle degrees of freedom. The latter are described approximately in terms of 
effective masses that depend on both the order parameter and the temperature, but 
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are independent of the symmetry breaking H term; these were presented in fig. |l|. 

The partition function then takes on a corresponding separable form and, as a 
consequence, it is possible to develop a simple method for performing a statistical 
sampling of the thermal equilibrium field configurations, including their time deriva- 
tives, at any temperature. This result, having a clear physical basis, is expected to 
be directly useful as practical means for initializing the dynamical simulations of the 
chiral field of the type carried out recently by several groups [^, ^, |TT], |T^, [I^, |TB[, 
thus making it easier to interpret the numerical results. 

Moreover, our specific illustrations provide useful insight into the equilibrium 
properties. In particular, it appears to be unrealistic to start the order parameter 
off with a value equal to zero. Indeed, at the critical temperature the most probable 
order parameter is much closer to its vacuum value f-,^ than to zero. The relation- 
ship between temperature and order parameter was summarized in fig. ^ and a more 
global impression of the distribution of the the order parameter (including its degree 
of misalignment) can be gained from the contour plots in fig. |[ 

Since the order parameter is thus very unlikely to vanish, the effective quasi- 
particle masses remain finite. Consequently, the statistical equilibrium distribution 
remains well behaved at all temperatures and the change from the "restored" phase 
to the normal one is fairly gradual. However, the finite size generally reduces the 
effective masses, thereby bringing the system closer to criticality. 

Of course, the statistical properties are of most practical interest at the relatively 
high temperatures characteristic of the initial stage of the high-energy collision. Once 
the chiral field has been initialized accordingly, any instabilities and associated ampli- 
fications will be automatically included in the dynamical propagation and the system 
can generally be expected to quickly move out of equilibrium. The equilibrium results 
can then provide a meaningful reference against which to analyze the deviation from 
equilibrium at any stage in the dynamical relaxation process. 

Additionally, we illustrated briefiy the equilibrium form of the correlation function 
which is an object of primary interest. Indeed, it is the correlation length that properly 
expresses the "domain size" governing the conjectured anomalous pion radiation. 
Essentially, what one would expect to see at the end is a stretched version of the 
initial correlation function since the long wave lengths are the most unstable and 
so will contribute in an ever larger propertion. This underscores the importance of 
starting out with chiral fields having physically reasonable correlation properties. To 
illustrate the use of the correlation function, we derived the rate at which real pion 
mesons are created by th field, and subsequently we calculated the dependence of the 
free pion strength on the temperature of the system. 

Finally, we sought to assess the degree of validity of our approach by subjecting 
sampled field configurations to the exact time evolution. This convenient means of 
testing suggested that the approximate treatment is of sufficient accuracy to be of 
practical utility. We therefore anticipate that it may find use in simulation studies, 
such as those exploring disoriented chiral condensates in high-energy collisions. 

The author is pleased to acknowledge helpful discussions with many colleagues, 
including P. Bedaque, L.P. Csernai, S. Gavin, J.I. Kapusta, Y. Kluger, V. Koch, L. 



22 



McLerran, R.D. Pisarski, R. Vogt, and X.N. Wang. 



23 



References 

[1] K. Rajagopal, in Quark-Gluon Plasma 2, ed. R. Hwa, World Scientific (1995) 

[2] S. Gavin, Nucl. Phys. A590 (1995) 163c 

[3] A.A. Anselm, Phys. Lett. B217 (1989) 169; 

A.A. Anselm and M.G. Ryskin, Phys. Lett. B266 (1991) 482 

[4] J.-P. Blaizot and A. Krzywicki, Phys. Rev. D46 (1992) 246 

[5] K. Rajagopal and F. Wilczek, Nucl. Phys. B399 (1993) 395 

[6] C.M.G. Lattes, Y. Fujimoto, and S. Hasegawa, Phys. Rep. 65 (1980) 151; 
J.R. Ren et al, Phys. Rev. D38 (1988) 1417; 
L.T. Baradzei et a/., Nucl. Phys. B370 (1992) 365 

[7] K.L. Kowalski and C.C. Taylor, Report CWRUTH-92-6 (1992); 

J.D. Bjorken, K.L. Kowalski, and C.C. Taylor, Report SLAC-PUB-6109 (1993) 

[8] K. Rajagopal and F. Wilczek, Nucl. Phys. B404 (1993) 577 

[9] S. Gavin, A. Gocksch, and R.D. Pisarski, Phys. Rev. Lett. 72 (1994) 2143 

[10] S. Gavin and B. Miiller, Phys. Lett. B329 (1994) 486 

[11] F. Cooper, Y. Kluger, E. Mottola, and J.P. Paz, Phys. Rev. D51 (1995) 2377; 
Y. Khigcr, F. Cooper, E. Mottola, J.P. Paz, and A. Kovner, Nucl. Phys. A590 
(1995) 581c 

[12] J.-P. Blaizot and A. Krzywicki, Phys. Rev. D50 (1994) 442 

[13] Z. Huang and X.-N. Wang, Phys. Rev. D49 (1994) 4335 

[14] M. Asakawa, Z. Huang, and X.-N. Wang, Phys. Rev. Lett. 74 (1995) 3126 

[15] D. Boyanovsky, H.J. de Vega, and R. Holman, Phys. Rev. D51 (1995) 734 

[16] L.P. Cscrnai and I.N. Mishustin, Phys. Rev. Lett. 74 (1995) 5005 

[17] S. Mrowczyhski and B. Miiller, Phys. Lett. 363B (1995) 1 

[18] M. Gell-Mann and M. Levy, Nuovo Cimento 16 (1960) 705 

[19] M. Aguilar-Benitez et a/., Phys. Rev. D50 (1994) 1173 

[20] G. Baym and G. Grinstein, Phys. Rev. D15 (1977) 2897 

[21] C. Itzykson and J.B. Zuber, Quantum Field Theory, McGraw-Hill, New York 
(1980) 

[22] E. Medeiros de Lima and J. Randrup, Nucl. Phys. A529 (1991) 115 

24 



[23] M. Toda, R. Kubo, and N. Saito, Statistical Physics I, Springer, New York (1992) 
[24] W.A. Bardccu and M. Moshc, Phys. Rev. D34 (1986) 1229 
[25] R.D. Pisarski, Phys. Rev. D52 (1995) R3773 

[26] J.I. Kapusta, Finite-temperature field theory, Cambridge University Press, New 
York (1989) 

[27] L.D. Landau and E.M. Lifshitz, Statistical Physics, Pergamon, New York (1980) 
[28] D. Horn and R. Silver, Ann. Phys. 66 (1971) 509 

[29] M. Gyulassy, S.K. Kauffmann, and L.W. Wilson, Phys. Rev. C20 (1979) 2267 

[30] R.D. Amado, F. Cannata, J. P. Dedonder, M.P. Locher, and Bin Shao, Phys. 
Rev. Lett. 72 (1994) 970 

[31] J. Bernstein and S. Dodelson, Phys. Rev. Lett. 66 (1991) 683 



25 




30 60 90 120 150 



Order parameter (|)o (MeV) 



Figure 1: Effective masses. 
The effective masses /xy (solid) and (dashed), as functions of the magnitude of the 
order parameter, 0o, for a range of temperatures: T — 0, 40, 80, 100, 122.63 {—Tc), 
160, 200, 240, 300, 400, 500 MeV, calculated in the thermodynamic limit where 
the box size is large, L oo. For a temperature above Tc, the two effective-mass 
curves start out at 0o = with degenerate values, whereas below they only exist 
if (po is sufficiently large. The corresponding starting points for /iy are connected by 
the dotted curve and, since nw is then nearly independent of T, only the curve for 
T — has been shown. The vertical arrow points to the vacuum value of the order 
parameter, </)vac = /tt = 92 MeV, and the free mass values = rrifj = 600 MeV/c^ 
and fi±^ = nij^ = 138 MeV/c^ arc indicated by the horizontal arrows. The locations 
of the corresponding points in the diagram are shown by the two solid symbols. 
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Figure 2: Effect of finite size on tlie effective masses. 
The effective masses are sliown as functions of tlie temperature for eitlier 0o = 0, wlien 
tlie 0(4) symmetry is restored, or for 0o = fn, tlie pliysical vacuum value. Tliree 
different scenarios fiave been considered: eitfier tfie tfiermo dynamic limit of large 
volume in which the quasi-particle spectrum is continuos (solid curve), corresponding 
to the scenario of fig. or a finite cubic box with side length L = 8 fm (dashed) or 
L = 5 f m (dotted), where the quasi-particle modes are quantized. 
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Figure 3: Fluctuations. 
The typical magnitude of the spatial fluctuations in the field strength, -< 5(j)^ 
as a function of the magnitude of the order parameter, for specified temperatures 
T: 122.63 (=Tc), 160, 200, 240 MeV. The system is enclosed in a cubic box of side 
length L ^ oo (solid curves), L = 8 fm (dashed), or L = 5 fm (dotted). The lower 
four curves display the fluctuations along the chiral direction of the order parameter 
</>, while the upper curves represent the total fluctuation. 
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Figure 4: Effective potential (part a: H = 0). 
The effective potential energy density Vt along the a axis for either the O (4) symmet- 
ric model having H — (part a) or with the adopted positive value of H (part b). The 
solid curves show the results for a number of temperatures: 100, 122.63 (=Tc), 160, 
200, 240 MeV. For T < Tc the effective potential curve starts at a certain minimum 
value of 00 between and v, These starting points are connected by the dotted curve, 
while the dashed curve shows the bare potential Vq obtained when fluctuations are 
neglected. The arrows point to the minima of Vt] for if > there is only a single min- 
imum (located at (Jq = fn), while for H = the degenerate ground-state minima form 
the surface of the 4-sphcre determined hj (f) = v. The value of the effective potential 
corresponding to other orientations of the order parameter can be easily obtained by 
noting that the directional dependence of Vt is given by — H<pQ cos Xo/T^^'c'. 
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Figure 4: Effective potential (part b: H > 0). 
The effective potential energy density Vt along the a axis for either the O (4) symmet- 
ric model having H — (part a) or with the adopted positive value of H (part b). The 
solid curves show the results for a number of temperatures: 100, 122.63 (=Tc), 160, 
200, 240 MeV. For T < Tc the effective potential curve starts at a certain minimum 
value of 00 between and v, These starting points are connected by the dotted curve, 
while the dashed curve shows the bare potential Vq obtained when fluctuations are 
neglected. The arrows point to the minima of Vt] for if > there is only a single min- 
imum (located at (Jq = fn), while for H = the degenerate ground-state minima form 
the surface of the 4-sphcre determined hj (f) = v. The value of the effective potential 
corresponding to other orientations of the order parameter can be easily obtained by 
noting that the directional dependence of Vt is given by — H<pQ cos Xo/T^^'c'. 
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Figure 5: Free energy (part a: H = 0). 
The free energy density Ft{^) along the a axis for either the 0(4) symmetric model 
having H = (part a) or with the adopted positive value of H (part b). The solid 
curves show the results for a number of temperatures. For T < Tc the curve starts at 
a certain minimum value of the magnitude 0o and these starting points are connected 
by the dotted curve, while the dashed curve shows the result obtained when the 
temperature is neglected. The arrows point to the minima of Ft; for if > there is 
only a single minimum, while for H = the degenerate minima form the surface of a 
4-sphere. For each temperature the location of the minima are indicated by the solid 
dots. (For H = the shallow secondary minima at 0o = are shown by open dots.) 
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Figure 5: Free energy (part b: H > 0). 
The free energy density along the a axis for either the 0(4) symmetric model 

having H = (part a) or with the adopted positive value of H (part b). The solid 
curves show the results for a number of temperatures. For T < Tc the curve starts at 
a certain minimum value of the magnitude 0o and these starting points are connected 
by the dotted curve, while the dashed curve shows the result obtained when the 
temperature is neglected. The arrows point to the minima of Ft; for if > there is 
only a single minimum, while for H = the degenerate minima form the surface of a 
4-sphere. For each temperature the location of the minima are indicated by the solid 
dots. (For H = the shallow secondary minima at 0o = are shown by open dots.) 
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Figure 6: Temperature dependence of the order parameter. 
The most probable value of (pQ, the magnitude of the order parameter, in the standard 
model where H > 0. In the limit of a very large box, L — oo (solid curve) 0o is 
constrained to the value for which the free energy density has its minimum (see fig. 
^) . The bars show the full width at half maximum of the thermal distribution of 0o 
in the system with L = 5 fm; those for L = 8 fm are about half that size (see fig. 0). 
The open dots show the centroids for H = for the box with L = 5 fm, while the 
light curve shows the corresponding result for a large box. 
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Figure 7: Distribution of the magnitude of the order parameter. 
For a cubic box of side length L = 8 fm is shown the probabihty density for the mag- 
nitude of the order parameter, P(0o)j for a range of specified temperatures T (indi- 
cated), obtained either by scahng the continuum result (solid curves) or by quantizing 
the quasi-particle modes (dashed curves). 
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Figure 8: The joint distribution P(0O)Xo)- 
The projection of the probabihty density -P(</>) onto the variables 4>o (the magnitude 
of the order parameter) and xo (the disahgnment angle) is displayed as a function 
of (Jo = 0ocosxo a-nd ttq = (/)osinxo) for a cubic box of side length L — 8 fm. 
For each temperature T, the solid dot indicates the location of the maximum of 
P{4>o, Xo) and the solid curve traces out the half-maximum contour, obtained by 
scaling the continuum results. For the temperatures 200 and 240 meV is indicated 
the corresponding result of quantizing the quasi-particle modes (dashed contours and 
open centroid dots). 
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Figure 9: Most probable effective masses. 
The most probable effective masses yUy (solid) and (dashed) are shown as functions 
of the temperature T, with the error bars indicating the full width at half maximum 
of the distribution resulting for a box with a side length of L = 8 fm. The results 
obtained by scaling the continuum results are indistinguishable from those of the 
corresponding quantized treatment. 
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Figure 10: Correlation function ai T — Tc. 
The reduced correlation function Cu S(f){ri) ■ 54>{r2) >- / ^ S(f)^ >- calculated 
at the critical temperature for various magnitudes of the order parameter, ranging 
from 00 = to the vacuum value, 0vac = /tt = 92 MeV. The most probable value at 
T = Tc is 00 f« 87 MeV. 



37 




1 2 3 4 5 



Separation r^2 (fnn) 



Figure 11: Temperature dependence of correlation function. 
The reduced correlation function C12 for a range of temperatures T, employing for 
each one the most probable value of (po, the magnitude of the order parameter. 
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Figure 12: Correlation length. 
The correlation length (the full width of C12 at half maximum) as a function of 
temperature, employing the most probable value of 4>o for each T. 
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Figure 13: Pion production rate. 
The strength function Pi_{m^c^) determining the overall rate at which vr'' mesons are 
being produced (see eq. (|55D) in a source in thermal equihbrium, as a function of 
its temperature T. This result is obtained by performing a Fourier transform of the 
quasi-particle correlation function associated with a given order parameter 0o and 
subsequently averaging over its distribution, Po(0o)- The overall normalization is 
arbitraty. 
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Figure 14: Early dynamics of the order parameter. 
The early trajectories of the order parameter is shown for a sample of eight config- 
urations, considering a box with L = 8 fm with a temperature of either 200 or 240 
MeV. The display is similar to that in fig. |] and the half-density contours (dashed) as 
well as the centroids (open dots) are those already given there. For each individual 
trajectory, the initial location is indicated by the sohd dot and the attached solid 
curve traces out the dynamical path up to the time t = 1 fm/c. 
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Figure 15: Time-averaged distribution of the order parameter. 
This figure illustrates the infiuence of the time evolution on the distribution of the 
magnitude of the order parameter, 0o, for a bow with L = 8 fm and for the tem- 
peratures of 180 and 240 MeV. The solid curves show the initial distribution of (pQ, 
as given by approximate statistical distribution ^^(^o) (see eq. (1^)). A sample of 
40 systems are then followed up to the time t = 10 fm/c and the order parameter 
is binned at regular time intervals throughout the evolution, leading to the dashed 
curves. 
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Figure 16: Time evolution of the correlation function. 
The pion correlation function 6*^2(512) at the temperature T = 240 McV. The dotted 
curve is the continuum limit [L 00) and the solid curve is the corresponding thermal 
result for a quantized finit box with a side length of L = 5 fm. The correlation function 
for a sample of ten initial configurations are shown by the short-dashed curves, and 
the long-dashed curves show the corresponding result after they have been propagated 
self-consistently up to the time t = 10 fm/c. The dashed curves have been obtained in 
two different ways: The curves that go up again result from aligning the separation r\2 
along one of the cartesian directions, while the other two are obtained for separations 
directed along a diagonal. The aligned curves have a periodicity equal to L, whereas 
the periodicity of the diagonal curves is -\/3 times larger. 
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